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Abstract — Cognitive  radio  can  dynamically  adapt  to  the 
available  spectrum  in  the  wireless  network.  Scheduling  and 
spectrum  allocation  are  tasks  affecting  the  performance  of 
cognitive  radio  wireless  network.  In  [1],  an  iterative  approach 
was  proposed  to  efficiently  compute  the  optimal  scheduling 
for  wireless  mesh  networks  with  single  channel  and  single 
radio.  The  optimal  scheduling  problem  is  decomposed  to 
a  sequence  of  small  optimization  problems  and  maximum 
weighted  independent  set  (MWIS)  problems,  and  both  of  them 
can  be  computed  quickly  even  for  large  networks.  For  example, 
the  optimal  scheduling  can  be  computed  for  the  mesh  network 
with  more  than  2,000  links  in  less  than  one  hour.  Here,  the 
iterative  algorithm  is  extended  to  the  cognitive  radio  wireless 
network  with  multi-channel  and  multi-radio.  Allowing  the 
schedule  problem  over  multi-channel  multi-radio  results  in 
higher  dimension  optimization  problem.  However,  the  proposed 
algorithm  can  obtain  the  optimal  spectrum  allocation  and  the 
schedule  quickly  for  moderate  size  of  networks.  Numerical 
experiments  show  that  the  optimal  throughput  is  achieved  when 
the  number  of  channels  is  one  or  two  more  than  the  number 
of  interfaces. 

I.  Introduction 

The  broadcast  nature  of  wireless  radio  causes  interference 
among  the  nodes  sharing  the  common  communication  me¬ 
dia.  To  reduce  the  interference,  fixed  spectrum  assignment 
policy  regulated  by  Federal  Communications  Commission 
(FCC)  assigns  a  fixed  portion  of  the  spectrum  to  a  license 
holder  or  a  wireless  service  for  exclusive  use  on  a  long  term 
basis.  However,  a  large  portion  of  the  assigned  spectrum  is 
sporadically  used  or  even  not  utilized  at  all,  which  leads  to 
the  scarcity  of  available  spectrum.  On  the  other  hand,  certain 
portions  of  the  spectrum  are  heavily  used,  such  as  2.4GHz 
band  and  5.0GHz  band. 

Dynamic  spectrum  allocation  (DSA)  is  a  good  way  to 
share  the  spectrum  among  users.  The  appearance  of  cogni¬ 
tive  radio  makes  DSA  feasible  where  cognitive  radios  can 
be  programmed  to  transmit  and  receive  on  a  variety  of  fre¬ 
quencies  and  use  different  transmission  access  technologies. 
However,  how  to  efficiently  share  the  available  spectrum  is 
a  challenge  in  the  cognitive  radio  wireless  networks  with 
multi-channel  and  multi-radio. 

* Visiting  scientist  under  a  National  Research  Council  Research  Associ- 
ateship  Award  at  the  Air  Force  Research  Laboratory. 


It  is  well  known  that  the  optimal  scheduling  problem  in 
wireless  network  is  NP-hard  even  when  single  channel  single 
radio  is  considered.  The  optimization  space  grows  exponen¬ 
tially  with  the  number  of  links.  For  example,  if  single  chan¬ 
nel  is  used,  then  one  must  maximize  over  a  polytope  with 
2l  extreme  points,  where  L  is  the  number  of  links.  More 
specifically,  we  define  an  assignment  to  be  a  specification  of 
which  links  are  transmitting  and  which  links  are  not  trans¬ 
mitting.  A  schedule  is  the  convex  sum  of  assignments.  The 
main  challenge  to  the  throughput  maximization  problem  has 
been  that  the  space  of  all  assignments  is  too  large.  In  order 
to  solve  the  scheduling  problem,  the  space  of  assignments 
must  be  reduced.  In  [1],  the  authors  have  developed  the 
iterative  techniques  to  compute  the  optimal  scheduling  for 
the  wireless  mesh  network  over  single  channel  even  when 
co-channel  interference  is  present.  The  idea  is  as  follows. 
Given  an  initial  set  of  randomly  selected  assignments  V, 
solve  the  throughput  optimization  problem  over  the  set  V 
and  obtain  the  Lagrange  multiplier  for  each  link.  Then,  find 
a  better  assignment  v+  by  solving  a  maximum  weighted 
independent  set  (MWIS)  problem.  If  v+  exists,  add  it  to 
the  set  V  and  go  back  to  solve  the  optimization  problem 
again.  Otherwise,  the  optimal  scheduling  is  obtained.  The 
optimality  and  convergence  property  are  proved  in  [1]. 

Allowing  the  schedule  optimization  to  optimize  over 
multi-channel  multi-radio  results  in  a  higher  dimension 
problem  than  when  the  optimization  is  over  only  single 
channel  single  radio.  However,  the  main  results  and  the 
algorithm  for  optimal  scheduling  from  [1]  are  still  available 
even  if  multi-channel  and  multi-radio  are  considered.  After 
solving  the  scheduling  problem,  the  assignment  specifies 
which  links  are  transmitting  and  their  transmission  channel. 
Therefore,  DSA  problem  is  also  solved  as  a  by-product.  In 
this  scheduling  problem,  the  dimension  of  MWIS  problem 
to  finding  a  new  assignment  significantly  increases  with  the 
number  of  channels.  Furthermore,  the  constraint  of  multi¬ 
radio  for  each  node  makes  the  MWIS  problem  even  harder. 
The  computation  complexity  of  solving  MWIS  determines 
the  complexity  of  the  scheduling  problem  over  multi-channel 
multi-radio. 

The  remainder  of  this  work  is  organized  as  follows. 
Section  II  briefly  introduces  some  related  work.  In  Sec- 


1 


tion  III,  notation  and  problem  definition  are  given.  Section 

IV  presents  the  iterative  approach  to  finding  the  spectrum 
allocation  and  optimal  scheduling  of  the  cognitive  radio 
networks  when  multi-channel  multi-radio  is  supported.  The 
approach  focuses  on  reducing  the  size  of  the  space  of 
assignments  and  finding  a  new  assignment.  Then,  Section 

V  shows  the  results  from  numerical  experiments.  Finally, 
concluding  remarks  are  given  in  Section  VI. 

II.  Related  Work 

Utility  optimization  of  wired  networks  was  pioneered  by 
the  seminal  works  of  Kelly  [2],  Low  [3].  The  same  frame¬ 
work  has  been  extended  to  ad  hoc  networks  in  [4],  [5].  The 
main  obstacle  facing  throughput  maximization  is  that  the 
size  over  which  the  optimization  is  performed  is  exponential 
with  the  number  of  links  [6]— [8]  even  when  single  channel 
single  radio  is  considered.  To  reduce  the  optimization  space, 
one  approach  is  to  arbitrarily  restrict  the  focus  to  a  subset 
of  all  assignments,  which  is  followed  in  [8].  However,  this 
approach  reduces  the  resulting  throughput.  Another  approach 
is  to  develop  a  heuristic  to  determine  the  subset  of  considered 
assignments  [9],  [10].  The  performance  of  these  methods 
typically  decreases  as  the  number  of  links  increases  [9],  [10]. 
Finally,  one  can  take  a  brute  force  approach,  and  consider 
the  entire  set  of  assignments.  This  approach  is  taken  in  [6], 
where,  due  to  computational  difficulties,  the  largest  network 
that  could  be  solved  has  only  15  links.  However,  the  authors 
in  [1]  has  developed  the  iterative  techniques  to  compute  the 
optimal  scheduling  for  the  wireless  mesh  network  with  more 
than  2,000  links. 

When  cognitive  radio  wireless  network  with  multi-channel 
single-radio  is  considered,  [11]  proposed  the  concept  of 
time-spectrum  block  and  the  protocols  based  on  heuristic 
method  to  allocate  these  blocks.  [12]  also  reduced  the 
optimization  space  by  constructing  the  subset  of  the  good 
assignments,  and  obtained  the  suboptimal  from  the  given 
assignments.  The  network  topology  considered  is  around  10 
links. 

Throughput  optimization  has  been  studied  for  wireless 
mesh  network  with  multi-channel  multi-radio.  [  1 3]— [  1 5] 
provided  the  theoretic  results  along  with  the  centralized 
algorithm.  Both  static  and  dynamic  spectrum  allocations  are 
investigated  in  these  works.  Due  to  the  large  size  of  opti¬ 
mization  space,  the  approximation  approaches  were  widely 
used  in  these  works,  and  the  proposed  algorithms  provided 
the  bounds  or  the  ratio  of  the  capacity  instead  of  the  optimal 
solution.  This  work  proposes  an  efficient  approach  to  finding 
the  spectrum  allocation  and  scheduling  for  cognitive  radio 
networks  with  multi-channel  multi-radio  by  extending  the 
work  in  [1]. 

III.  Notation  and  Problem  Definition 

A  router-to-router  flow  is  denoted  by  <j>,  with  $  denoting 
the  set  of  all  such  flows.  To  improve  presentation,  it  is 


assumed  that  all  flows  use  a  single  path,  however,  the 
extension  to  multipath  is  straightforward.  The  data  rate  of 
flow  </>  is  denoted  by  /^,  and  the  path  followed  by  flow 
</>  is  denoted  by  P  (cf).  The  set  of  all  considered  paths  is 
V .  Using  this  notation,  the  total  data  rate  sent  over  link 
£  is  E{0| where  i'/’l21  €  P{4>)}  is  the  set  of 
flows  that  cross  link  x.  All  links  are  directed;  bidirectional 
communication  between  two  nodes  is  represented  by  two 
directed  links. 

Define  an  assignment  to  be  a  vector  v  = 

[  V\  •••  Vl  ],  where  there  are  L  links  in  the  network 
and  where  vx  £  {0, 1}  with  vx  =  1  implying  that 
link  x  is  transmitting  during  assignment  v.  The  set  of 
considered  assignments  is  denoted  by  V,  while  the  set  of 
all  assignments  is  denoted  by  V.  In  this  single  channel 
case  where  vx  £  {0, 1},  V  has  2L  assignments.  The  size  of 

V  is  the  main  challenge  facing  optimal  scheduling.  Thus, 
typically,  we  only  consider  a  subset  of  all  assignments,  i.e., 

V  cy. 

Suppose  that  the  available  bandwidth  B  is  divided  evenly 
into  K  orthogonal  channels,  each  having  bandwidth  A B  = 
B/K,  sufficiently  small  to  be  accurately  approximately  by 
a  flat  fading  coefficient.  For  simplicity,  assume  all  channels 
have  the  same  fading  coefficient.  Also,  assume  that  there  are 
L  links  in  the  cognitive  radio  wireless  network,  and  each 
node  x  has  A'IX  cognitive  radio  interfaces.  For  a  given  node 
rid,  all  links  incident  to  node  rid  are  denoted  by  E(rid). 

When  considering  multiple  channels,  the  concept  of  a 
logical  link  arises.  For  example,  a  physical  link  may  transmit 
at  K  orthogonal  channels.  This  can  be  modeled  as  K  logical 
links,  each  with  a  different  channel,  but  between  the  same 
pair  of  physical  nodes.  In  this  case,  the  logical  link  Xk 
would  refer  to  the  physical  link  x  transmitting  at  the  A'th 
channel.  Then,  an  assignment  specifies  which  logical  links 
are  transmitting.  We  drop  the  terms  physical  and  logical 
unless  it  is  unclear  from  context. 

A  schedule  is  defined  as  a  convex  combination  of  assign¬ 
ments.  For  example,  in  the  simple  case  where  all  links  can 
only  transmit  at  a  single  channel,  then  an  assignment  can 
be  written  as  v  €  {0, \}L ,  where  there  are  L  links  in  the 
network.  In  this  case,  vx  =  1  implies  that  link  x  is  trans¬ 
mitting  during  assignment  v.  A  schedule  can  be  represented 
by  a  set  {av}  where  av  >  0  and  av  <  1-  Letting 
R  (v,  x)  be  the  data  rate  over  link  x  during  assignment  v,  the 
average  data  rate  over  link  x  provided  by  schedule  {av}  is 
avR  (v,  x).  More  details  on  the  single  channel  case  can 
be  found  in  [1].  Since  links  are  permitted  to  use  K  channels; 
hence  an  assignment  can  be  written  as  v  £  {0,  l}LxK , 
where  vXk  =  1  implies  that  the  physical  link  x  is  transmitting 
at  channel  k.  Due  to  the  large  size  of  {0,  l}LxA,  we  consider 
subsets  V  C  {0,  l}LxA,  where  V  is  referred  to  as  the  set  of 
considered  assignments. 

The  data  rate  across  logical  link  Xk  during  assignment  v 
is  denoted  by  R{y,Xk).  A  simple  binary  approximation  is 
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used  to  define  R(y,Xk )•  Specifically, 


{RXk  if  vVj  =  0  for  all  y:j 
Vi  7 k  xk,vXk  =  1, 

0  otherwise, 

(1) 

where  x(xk)  is  the  set  of  logical  links  that  are  in  conflict 
with  Xk,  i.e.,  yj  €  x  ixk)  if  simultaneous  transmissions 
over  logical  link  Xk  and  logical  link  yj  are  not  possible. 
RXk  is  the  nominal  data  rate  over  logical  link  Xk .  Note 
that  this  definition  of  R{v,Xk)  neglects  the  possibility  that 
the  aggregate  interference  from  transmissions  over  several 
links  not  in  \  (x k )  can  result  in  a  transmission  failure  over 
logical  link  xk.  A  technique  to  account  for  such  interference 
is  discussed  in  Section  IV-B.  All  computations  in  this  work 
will  use  this  technique,  and  hence  the  computed  throughputs 
account  for  multiple  interferers. 

The  set  of  conflicting  links,  %  ( Xk ),  depends  on  the 
communication  model.  Arguably,  the  SINR '.binary  commu¬ 
nication  model  is  the  most  relevant  and  is  the  model  used 
in  this  work.  Here,  the  SINR  is  defined  as 


SINR  (xk,yj) 


Hxk,xkPxk 
Hy^XkPyj  +  A/” 


where  HXh  Xk  is  the  channel  gain  across  link  xk,  Hy, 
is  the  channel  gain  from  the  transmitter  of  link  yj  to  the 
receiver  of  link  xk,  pXk  and  pyj  are  the  transmission  powers, 
and  A f  is  the  channel  noise.  Let  us  define  T(xk)  be  the  SINR 
threshold  to  achieve  the  desired  data  rate  RXk.  Under  the 
SINR  binary  communication  model,  the  set  of  conflicting 
links,  x(xk)  is  defined  by 


X(xk)  =  {Vj  \SINR(xk,yj)  <  r(xfe)  (2) 
or  SINR  (yj,Xk)  <  T (yj)}. 


Suppose  L  links  share  K  channels  and  each  node  x 
has  Mx  cognitive  radio  interfaces.  With  this  notation,  the 
schedule  optimization  problem  can  be  described  as 


maxG  (f) 

C*,f 

subject  to: 


(3a) 


K 


E  u<Y,  av  R  (v,  Xk)  for  each  link 

{<t>\x€P(<t>)}  v£V  k= 1 


E«v<i 

vev 

0  <  av  for  each  v  €  V. 


(3b) 

(3c) 

(3d) 


where  f  is  the  vector  of  flow  rates.  The  function  G  is  referred 
to  as  the  throughput  metric.  Several  different  throughput 
metrics  are  possible.  In  some  cases,  the  throughput  metric 
is  the  network  utility  G(f)  =  where 

is  the  utility  function  for  flow  </>.  Popular  utility  func¬ 
tions  include  U<j,{f)  =  u^log(/)  [2],  [16],  [17]  and 


Algorithm  1  Computing  an  Optimal  Schedule 
1:  Select  an  initial  set  of  assignments  V(0),  set  i  =  0. 

2:  Solve  (3)  for  V  =  V(z)  and  compute  fi(i)  and  A (z), 
the  Lagrange  multipliers  associated  with  constraints  (3b) 
and  (3c),  respectively. 

3:  Search  for  an  assignment  v+  that  solves 

L  K 

v+  €  argmax  (4) 

x—1  k— 1 

4:  if  the  assignment  v+  €  V(z)  or 
Ex=i  Vx(i)  Ek= iR(v+,xk)  <  A (i)  then 

5:  stop,  the  optimal  schedule  has  been  found. 

6:  else 

7:  set  V  (z  +  1)  =  V  (z)  U  v+,  set  z  =  z  +  1,  and  go  to 

Step  2. 

8:  end  if 


Ucf>  (/)  =  w^f1  1 1  (1  —  7)  [18],  where  w $  is  the  admin¬ 
istrative  weight.  Another  widely  used  throughput  metric  is 
G  (f)  =  [8],  If  G(f)  =  min ^  then 

(3)  can  be  written  as  a  linear  programming  problem,  which 
is  extensively  studied  in  [19].  This  work  specifically,  focuses 
on  the  case  when  G  (f)  =  rnin^g^, 

IV.  Optimal  Scheduling 

In  theory,  (3)  is  solvable.  However,  there  is  a  significant 
computational  challenge  in  that  if  V  is  the  set  of  all  assign¬ 
ments,  then  the  vector  a.  has  2LK  elements.  Thus,  the  size 
of  the  space  over  which  the  optimization  is  performed  must 
be  reduced.  This  idea  of  considering  a  reduced  space  was 
considered  in  [8]  and  [9],  however,  the  space  was  constructed 
arbitrarily.  In  this  work  the  space  is  constructed  so  that  the 
throughput  found  by  optimizing  over  the  reduced  space  is 
the  same  throughput  found  by  optimizing  over  the  entire 
space. 

A.  Finding  New  Assignments 

In  order  to  find  a  new  assignment,  (4)  must  be  solved.  It 
is  well  known  that  solving  (4)  is  equivalent  to  solving  the 
graph  theoretic  problem  known  as  the  maximum  weighted 
independent  set  (MWIS)  problem  (e.g.,  see  [1]).  In  the 
worst  case,  the  MWIS  problem  is  NP-hard.  However,  this 
does  not  mean  that  the  MWIS  problem  is  computational 
intractable  in  all  cases.  In  fact,  [20]  showed  that  (4)  can  be 
solved  efficiently.  For  example,  [20]  found  that  in  practical 
networks,  (4)  can  be  solved  in  about  one  second  for  a 
topology  with  2048  links  where  single  channel  single  radio 
is  considered. 

One  approach  to  solving  (4)  is  to  use  a  generic  integer 
linear  programming  (ILP)  solver.  The  MWIS  problem  can 
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be  formulated  as  an  ILP  via 

L  K 

max 

Xk  ^Efc  (5a) 

x=l  k=l 

subject  to:  vXk  +  vy.  <  1  if  yj  e  X  ( %k ) 

K 

E  E  vXk  <  M.y  for  each  node  y  (5b) 

x£E(y)  k—1 

Vxk  €  {0,  1}  . 

where  x  €  E(y)  denotes  the  links  incident  to  node  y 
and  (5b)  restricts  that  the  active  links  incident  to  node  y 
must  be  no  more  than  the  number  of  cognitive  radios  in 
node  y.  The  time  required  to  solve  this  problem  can  be 
significantly  decreased  if  a  greedy  clique  decomposition  of 
the  conflict  graph  is  performed.  Specifically,  a  set  of  cliques 
{Qi,  i  —  1,2,  ...Q}  are  found  such  that  if  y,  €  x  (xk ).  then 
there  is  a  clique  Qi  such  that  Xk  €  Qi  and  yj  €  Qi.  Then, 
(5a)  becomes 

L  K 

max  EX  E  R-'x  VXk  (6) 

x—1  k—1 

subject  to:  E^  vXk  <  1  for  i  =  1, 2, Q 

xk^Qi 

K 

E  E  vXk  <  My  for  each  node  y 

xeE(y)  k= 1 

VXk  S  {0, 1}  . 

B.  Removing  Multi-Conflict 

Binary  communication  models  have  the  drawback  that 
they  do  not  directly  account  for  the  impact  of  the  aggregate 
of  interference  from  several  simultaneously  transmitting 
links.  For  example,  suppose  that  links  y:j  ^  x  (xk),  Zi 
X(xk),  Vj  ^  xizi)>  thus,  according  to  the  binary  conflict 
model,  any  pair  of  these  three  links  can  simultaneously 
transmit.  Flowever,  it  is  possible  that  the  aggregate  of  the 
interference  of  link  y:j  and  link  zi  transmitting  is  strong 
enough  that  communication  across  link  Xk  is  not  possible. 
Flowever,  binary  conflict  models  do  not  account  for  this  type 
of  interference.  Consequently,  a  schedule  constructed  based 
a  binary  conflict  model  may  perform  poorly  when  applied 
in  the  physical  model. 

One  way  to  account  for  these  multi-conflict  is  as  follows. 
Let  v+  be  a  new  assignment  found  by  solving  (6).  Check 
whether  this  assignment  has  any  multi-conflict,  by  determin¬ 
ing  whether  there  is  a  Xk  with  u+  =  1  and 

_ Hxk,xkPxk _  ,  p/  \ 

x^k  n  \r  ^  1  \Xk>- 

Z^y=i,y9kx  iEj=i  nVj,xkPyjvyj  +-/v 

If  such  a  Xk  exists,  then  add  the  found  multi-conflict  into 
the  ILP  (6)  as  follows.  Let  C  =  =  l},  that  is,  C  is 

the  set  of  links  that  make  up  the  multi-conflict.  Intuitively, 


C  should  be  the  smallest  set  of  links  that  form  the  multi¬ 
conflict  at  link  Xk .  This  multi-conflict  can  be  removed  from 
consideration  by  using  the  following  ILP  problem  to  search 
for  new  assignments, 

L  K 

max  R-'x  VXk  (7a) 

x=l  k—1 

subject  to:  ^  vXk  <  1  for  i  =  1,  2, ...,  Q 

xktQi 

K 

E  E*.  <  My  for  each  node  y 

x£E(y)  k=  1 

E <  ici  - 1  (7b) 

Xk  GC- 

VXk  ^  {bi  1 

where  \C\  is  the  number  of  elements  in  C.  Again,  it  is 
determined  whether  the  assignment  found  by  solving  (7) 
has  multi-conflict  and  if  so,  another  constraint  such  as  (7b)  is 
added  to  the  ILP  This  process  continues  until  an  assignment 
is  found  that  does  not  have  a  multi-conflict. 

V.  Numerical  Experiments 
A.  Computational  Experimental  Set-up 

In  order  to  determine  the  performance  of  the  various  tech¬ 
niques  described  above,  this  investigation  employed  the  re¬ 
alistic  urban  mesh  networks  constructed  by  the  UDelModels 

[21] .  Along  with  a  realistic  propagation  simulator,  the  UDel 
Models  include  a  map  builder,  a  realistic  mobility  simulator, 
and  a  large  collection  of  data  and  trace  files.  The  propagation 
simulator  is  based  on  ray-tracing  and  accounts  for  reflections 
off  of  the  ground  and  off  of  buildings,  transmission  through 
building  walls,  and  diffraction  around  and  over  buildings 

[22] .  It  also  accounts  for  the  impact  that  different  materials 
have  on  reflections  off  of  walls  and  transmission  through 
walls.  Data  sets  for  several  urban  areas  are  available  online. 

For  this  investigation,  the  topologies  are  generated  as 
follows.  Each  topology  was  based  on  a  different  6x6  city 
block  region  that  was  randomly  located  within  a  2  km2 
region  of  downtown  Chicago.  Various  node  densities  were 
investigated.  Specifically,  the  number  of  gateways1  was  1,  3, 
and  6,  and  the  number  of  wireless  routers  was  18,  36,  54,  72, 
and  90.  The  wireless  routers  and  gateways  were  uniformly 
distributed  throughout  the  6x6  city  block  region. 

In  these  experiments,  all  traffic  flowed  from  the  gateways 
to  destinations  (i.e.,  downstream  traffic),  where  each  mesh 
router  in  the  topology  was  a  destination  of  a  flow.  The 
routing  was  a  single  path,  least  hop  routing,  where  each 
link  had  a  receiver  signal  strength  of  at  least  -75  dBm  when 
the  transmission  power  was  15  dBm. 

'Gateways  have  both  wired  and  wireless  interfaces,  whereas  wireless 
routers  only  have  wireless  interfaces. 
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Fig.  1.  Variation  in  the  computed  throughput  as  assignments  are  added 
where  the  throughput  is  min^e$  The  plot  is  for  a  topology  with  3  GW 
and  90  wireless  routers  where  the  number  of  channels  is  7  and  the  number 
of  radios  at  each  node  is  4. 

802.11a  data  rates  were  used.  Furthermore,  it  was  assume 
that  the  relationship  between  SNR  and  bit-rate  is  the  same  as 
the  relationship  between  SINR  and  bit-rate.  The  throughput 
metric  used  is  G  (f)  =  min^g.j  /^.  The  computations  below 
were  performed  on  a  2.4MHz  AMD  FX-53  processor  with 
8GB  RAM.  Programs  were  written  in  Matlab  and  used 
CPLEX  vlO  for  solving  LP  and  ILP  problems. 

B.  Number  of  Iterations  until  Algorithm  1  Stops 

Figure  1  shows  how,  in  the  topology  with  3  GW  and  90 
wireless  routers  where  the  number  of  channels  is  7  and  the 
number  of  radios  at  each  node  is  4,  the  throughput  increases 
as  the  more  assignments  are  added.  The  point  of  maximum 
throughput  occurs  when  the  stopping  condition  specified  in 
Algorithm  1  is  met.  Thus,  in  this  case,  Algorithm  1  stopped 
after  76  iterations  when  the  throughput  metric  was  G  (f)  = 
min^  fy. 

Note  that  only  one  assignment  is  added  at  each  iteration. 
Thus,  the  maximum  number  of  elements  in  V  is  the  number 
of  assignments  in  the  initial  set  of  assignments  plus  the 
number  of  iterations  required  by  Algorithm  1.  Hence,  we 
have  achieved  the  goal  of  determining  the  solution  to  (3)  for 
V  =  V  by  computing  the  solution  to  (3)  for  a  small  set  V. 
The  complexity  of  solving  linear  and  nonlinear  optimization 
problems  is  well  known,  and  is  not  investigated  here. 

C.  Multi-Channel  and  Multiple  Cognitive  Radios 

The  proposed  framework  easily  covers  the  multi-channel 
multi-radio  scenario.  Assume  that  all  nodes  have  the  same 
number  of  interfaces.  Figure  2  shows  the  relationships 
among  the  throughput,  the  number  of  channels  and  the 
number  of  interfaces.  The  y-axis  is  the  throughput  ratio  over 
the  optimal  throughput  with  single  channel  and  single  radio. 
These  plots  are  for  a  wireless  network  with  3  GWs  and  36 
wireless  routers.  A  few  comments  are  in  order. 

•  If  there  are  n  channels  available  and  n  radios  at  each 
node,  the  throughput  increases  with  a  factor  of  n 
because  all  channel  are  same  as  each  other. 


Number  of  Channels 


(a) 


(b) 


Fig.  2.  (a).  Throughput  versus  the  number  of  channels  with  different 

number  of  interfaces,  (b).  Throughput  versus  the  number  of  interfaces  with 
different  number  of  channels.  These  plots  are  for  a  wireless  network  with 
3  GWs  and  36  wireless  routers. 

•  The  optimal  throughput  is  achieved  when  the  number 
of  channels  is  one  or  two  more  than  the  number  radios. 
The  throughput  is  not  improved  any  more  even  if  the 
number  of  channel  increases  further. 

•  The  throughput  of  two  channels  and  single  radio  is 
better  than  that  of  single  channel  and  single  radio. 
Similar  results  are  correct  for  n  +  1  channels  and  n 
radios  case. 

The  results  are  also  true  for  other  generated  topologies. 

VI.  Conclusion 

In  this  paper,  we  presented  the  problem  of  joint  spectrum 
allocation  and  scheduling  in  multi-radio  multi-channel  cog¬ 
nitive  radio  wireless  networks.  The  iterative  algorithm  pro¬ 
posed  in  [1]  was  extended  to  solve  the  optimization  problem 
with  higher  dimension.  Specifically,  with  the  introduction 
of  logical  links,  the  dimension  of  the  MWIS  problem  is 
increased  by  I\  times  where  K  is  the  number  of  channels, 
while  the  optimization  problem  (3)  over  a  small  set  of 
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assignments  still  has  the  same  dimension.  With  this  iterative 
algorithm,  we  can  efficiently  compute  the  joint  spectrum 
allocation  and  scheduling  for  the  wireless  networks  with  a 
few  hundred  links.  In  addition,  numerical  experiments  show 
that  the  optimal  throughput  is  achieved  when  the  number  of 
channels  is  one  or  two  more  than  the  number  of  interfaces. 

Disclaimer 

This  work  was  performed  while  the  first  author  held  a 
National  Research  Council  Research  Associateship  Award 
at  the  Air  Force  Research  Laboratory. 

The  views  and  conclusions  contained  in  this  document 
are  those  of  the  authors  and  should  not  be  interpreted 
as  representing  the  official  policies,  either  expressed  or 
implied,  of  the  Air  Force  Research  Laboratory  or  the  U. 
S.  Government. 
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